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Abstract 



m 

We propose an intensity-dependent pion-nucleon coupling Hamiltonian within a 
\ unitary multiparticle-production model of the Auerbach-Avin-Blankenbecler-Sugar 

(AABS) type in which the pion field is represented by the thermal-density ma- 
trix. Using this Hamiltonian, we explain the appearance of the negative-binomial 
(NB) distribution for pions and the well-known empirical relation, the so-called 



Wroblewski relation, in which the dispersion D of the pion-multiplicity distribution 



is linearly related to the average multiplicity < n > : D = A < n > +B, with the 
coefficient A < 1. The Hamiltonian of our model is expressed linearly in terms of 
the generators of the SU(1, 1) group. We also find the generating function for the 
pion field, which reduces to the generating function of the NB distribution in the 
limit T -> 0. 
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1 Introduction 



During the last years a considerable amount of experimental information has been 
accumulated on multiplicity distributions of charged particles produced in pp and 
pp collisions in the centre- of-mass energy range from lOGeV to lSOOGeV . Mea- 
surements in the regime of several hundred GeV [1] have shown the violation of 
the Koba-Nielsen-Olesen (KNO) scaling [2], which was previously observed in the 
ISR cm. energy range from 11 to 63 GeV [3]. The violation of the KNO scaling is 
characterized by an enhancement of high-multiplicity events leading to a broadening 
of the multiplicity distribution with energy. 

The shape of the multiplicity distribution may be described by either its 
C moments, C q = {n q )/(n) q , or by its central moments (higher-order dispersions), 

D q = {{n— (n)) q ) 1 ^ q , q = 2, 3 The exact KNO scaling implies that all C q moments 

are energy independent. The C moments appear to be energy independent only at 
energies below lOOGeV. It can also been shown [4] that the KNO scaling leads to a 
generalized Wroblewski relation [5] 

D q = A q (n)-B q , (1) 

with the energy- independent coefficients A q and B q . The pp and pp inelastic data 
below lOOGeV also show a linear dependence of the dispersion on the average number 
of charged particles, but with the coefficients A q and B q that are approximately equal 
within errors. 

The fact that the dispersion of the multiplicity distribution grows linearly with 
(n) implies that the elementary Poisson distribution resulting from the independent 
emission of particles is ruled out. 

The total multiplicity distribution P n of charged particles for a wide range of 

energies (22 — QWGeV) is found to be well described by a negative-binomial (NB) 
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distribution [1,6] that belongs to a large class of compound Poisson distributions [7]. 
It is a two-step process [8] with two free parameters: the average number of charged 
particles (n) and the parameter k which affects the shape (width) of the distribution. 
The parameter k is also related to the dispersion D = D 2 by the relation 

so that the observed broadening of the normalized multiplicity distribution with 
increasing energy implies a decrease of the parameter k with energy. The KNO 
scaling requires k to be constant. 

Although the NB distribution gives information on the structure of correlation 
functions in multiparticle production, the question still remains whether its clan- 
structure interpretation is simply a new parametrization of the data or has deeper 
physical insight [9]. Measurements of multiplicity distributions in pp collisions at 
TeV energies [10] have recently shown that their shape is clearly different from that 
of the NB distribution. The distributions display the so-called medium-multiplicity 
"shoulder", with a shape qualitatively similar to that of the UA5 900GeV and UAl 
distributions [11]. A satisfactory explanation of this effect is still lacking [12]. 

In this paper we propose another approach to multiplicity distributions based 
on a unitary eikonal model with a pion-field thermal- density operator given in terms 
of an effective intensity-dependent pion-nucleon coupling Hamiltonian. We assume 
that the system of produced hadronic matter behaves as a hadron gas in thermody- 
namical equilibrium at the temperature T before the hadrons themselves decouple 
(freezing-out) and decay producing observable particles in the detector. 

The paper is organized as follows. In Sect. 2 we explain the basic ideas of our 
unitary eikonal model with a pion-field thermal-density operator. A discussion of 
the Wroblewski relation and the NB distribution are presented in Sect. 3. Finally, 



in Sect. 4 we draw conclusions and make remarks on a possible extension of the 
model to include two-pion correlations in the effective pion-nucleon Hamiltonian. 
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2 Description of the model 



At present accelerator energies the number of secondary particles (mostly pions) 
produced in hadron-hadron collisions is large enough, so that the statistical ap- 
proach to particle production becomes reasonable . Most of the properties of pions 
produced in high-energy hadron-hadron collisions can be expressed simply in terms 
of a pion-field density operator. We neglect difficulties associated with isospin and 
only consider the production of isoscalar "pions". In high-energy collisions, most of 
the pions are produced in the central region. In this region the energy - momentum 
conservation has a minor effect if the transverse momenta of the pions are limited 
by the dynamics. 

A. The AABS model 

In a class of unitary eikonal models (AABS models) [13] formulated a long 
time ago, the incident hadrons propagate through the interaction region, without 
making significant changes in their longitudinal momenta (leading-particle effect). 
Only the part W = Kyfs of the total cm. energy y/s in every concrete event is 
avaliable for particle production, where K is the inelasticity < K < 1 . 

In the AABS type of models, the scattering operator S is diagonal in the 
rapidity difference Y = ln(s/m 2 ) and in the relative impact parameter B of the two 
incident hadrons. The initial-state vector for the pion field is S(Y,B) | >, where 
the vacuum state | > for pions is in fact a state containing two incident hadrons. 

The n-pion production amplitude for n > 1 is given by 

tT n (Y,B;k 1 ...k n ) = 2s(k l ...k n | S(Y,B) | 0). (3) 

The square of the n-pion production amplitude can be written in the form 

| T n (Y, B;h... k n ) \ 2 = As 2 Tr{p(Y, B)\k 1 ... k n ) (h . . . k n | }, (4) 
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where the pion-density operator p(Y, B) is defined as 

p(Y,B) = S(Y,B)\0}(0\S\Y,B). (5) 

The square of the elastic scattering amplitude is then the matrix element of p(Y, B) 
between the states with no pions, i.e. (0 | p(Y,B) | 0). 

In terms of the pion-number operator 

N = a l a k = J2^k, k= (u k , k), (6) 

k k 

the square of the S'-matrix element when no pions are emitted can also be written 
in the form 

| (0 | S(Y, B) | 0) | 2 = Tr{p{Y, B) : e~* :} = e~ n(Y ^ . (7) 

Here : : indicates the operation of normal ordering and fl(Y, B) is the usual eikonal 
function (or the opacity function) of the geometrical model [14]. The connection 
with the inelastic cross section and the exclusive cross section for production of n 
pions is then 

a mel (Y, B) = 1 - e~ n ^ S \ (8) 

and for n > 1, it is 

a n (Y,B)=Tr{p(Y,B):—e- N :}. (9) 

In terms of a normalized pion-multiplicity distribution at each impact parameter, 
P n (Y,B) = <J n (Y, B) I <j ine [(Y, B) , the observed complete multiplicity distribution 

— * — * 

P n (Y) is obtained by summing P n (Y, B) over all impact parameters B with the 
weight function Q(Y, B) = a ine i(Y, B) / a inel (Y) , i.e. 

P n (Y) = J d 2 BQ{Y, B)P n {Y, B). (10) 

The first-order moment of P n (Y) gives the average multiplicity 

(n) = Y, nP n{y) = J d 2 BQ(Y, B)n(Y, B). (11) 
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The higher-order moments of P n (Y) give information on the dynamical fluctuations 
from (n) and also on the multiparticle correlations. All these higher-order moments 
can be obtained by differentiation from the pion-generating function 

G(z) = *£z n P n (Y) = I d 2 BQ(Y,B)G(Y,B;z), (12) 

where 

G(Y,B;z)=Tr{p(Y,B)z*} (13) 
is the pion-generating function in _B-space . Thus the normalized factorial moments 



F q are 



(n(n-l) („-,+ !)> = fOm 
q (n)i x ' dz« v 1 



and the normalized cumulant moments K n are 



^W-^. (15) 
These moments are related to each other by the formula 

^ = e(V)^-^- (16) 

For the Poisson distribution, all the normalized factorial moments are identically 
equal to unity and all cumulants vanish for q > 1. 

We shall be concern mostly with the q = 2 moments, which are directly related 
to the dispersion D: 

M + 1 . ( _^ + 1 __L (17) 

B. Thermal-density operator for the pion field 

— * 

The operator | 0)(0 | appearing in the definition of p(Y, B) represents the 
density operator p{vac) for the pion-field vacuum state. 
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The density operator for a pion field in thermal equilibrium at the temperature 

T is 
where 

H = £>*(4 a * + A )> ( 19 ) 

k 

InZ = -f3\J2^k-J2 ln ( 1 - e ~^ k )- 

k k 

The quantity A J2k ^k — (0 ] H \ 0)—J2k ^fc(0 | N k | 0) represents the lowest possible 
energy of the pion system in the leading - particle environment. The " zero-point 



energy " corresponds to A = \. If the energies Uk = yk 2 + m\ of the pion gas in 
the volume V are closely spaced, the summation over k is replaced by an integral: 
E k ^VJd 3 k/2u; k . 

Note that piyac) = pr=o- The mean number of thermal (chaotic) pions is 

= < 20 > 

k e 1 

= J2 n Tk- 

k 

Owing to the interaction of pions with the nucleon field the density operator 
Pt is transformed by means of the unitary ^-matrix into 

p T (Y,B) = S(Y,B)pr&(Y,B) (21) 

1 

Z 



_ e -I3H(Y,B) 



where 

H(Y, B) = S{Y, B)H &(Y, B) (22) 

is regarded as an effective Hamiltonian describing the pion system in the interaction 
with the leading - particle system. 



Now, taking into account an old observation of Golab-Meyer and Ruijgrok [15] 
that the Wroblewski relation can be satisfied for all energies if the square of the pion- 
nucleon coupling constant increases linearly with the mean number of pions (n), we 
propose the following form of the effective pion-nucleon coupling Hamiltonian: 

H(Y,B) = 5> fc (y, B)(N k + A) + g k (Y, B)(a k ^N k + 2A - 1 + h.c.)] (23) 
k 

= £if fc (Y,B), 

k 

where e k (Y, B) = u k +Ag k (Y, B). The interaction part of the Hamiltonian H k for the 
k mode is no longer linear in the pion-field variables a k and represents an intensity- 
dependent coupling [16]. It is also easy to see that the operators 

K (k) = N k + X, (24) 
K_(h) = a k y/N k + 2X-1, 
K+(k) = y/N k + 2X-lal 

form the standard Holstein-Primakoff [17] realizations of the su(l, 1) Lie algebra, 
the Casimir operator of which is 

C k = K 2 (k) - l -[K + {k)K_(k) + K_(k)K + (k)} = A(A - l)I k . (25) 

The Hamiltonian H k (Y, B) = H k is thus a linear combination of the generators of 
the SU(1, 1) group: 

H k = e k K (k) + g k [K + {k) + K_{k)]. (26) 

— * 

The corresponding S-matrix which diagonalizes the Hamiltonian H(Y, B) is 

S{Y,B) = \{S k {Y,B), (27) 

k 
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where 

S k (Y, B) = exp{-9 k (Y, B)[K + {k) - K_(k)}}, (28) 

with 

tM t (Y,B) = 2 -^§. (29) 
e k {Y,B) 

Since the dependence on the variables Y, B is contained only in the hyperbolic angle 
9 k (Y, B), we shall from now on assume this dependence whenever we write 9 k . 

It is easy to see that the initial-state vector for the pion field, S(Y, B) | 0), 
factorizes in the /c-space as 

S(Y,B) \0) = H(S k (Y,B) |0*», (30) 
k 

with 

S k (Y,B)\0 k ) = (1 -t^0 A E(-^) Wfc ( r(n frm^ 1/2 ' ^ ^ 

n h n fc !r(2A) 

= I 9k), 

where | n k ) = {n k \)~ l l 2 (a\) nk \ k ). In the same way we find that the pion thermal- 
density operator pt(Y,B) is also factorized as 

p T (Y,B) = l[p T (e k ), (32) 

k 



with 



pAOk) = ^j2 e ~ Mnk+x) I n k ,e k )(n k ,e k \, (33) 



z k 



where | n k ,9 k ) = S k (Y,B) \ n k ). The states | n k ,9 k ) form a complete orthonormal 
set of eigenvectors of the k-mode Hamiltonian H k , i.e. 

H k I n k ,9 k ) = u k (n k + A) | n k ,9 k ), (34) 
J2 I n k,0k)(n k ,9 k | = J, (35) 
(n k ,9 k | m k ,9 k ) = 5 n>m . (36) 
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3 Pion-generating function and its moments 

The average multiplicity n T (Y, B), the dispersion d^(Y, B) and all higher-order mo- 
ments at the temperature T in B space, 

4(F, B) = Tr{p T (Y, B)N«}, q = 1, 2, . . . (37) 

can be obtained by differentiation from the pion-generating function 

G T (Y,B;z) = l[G T (6 k ;z), (38) 
k 

where 

G T (9 k ;z)=Tr{p T (9 k )z^}. (39) 

After performing a certain amount of straightforward algebraic manipulations, 
we find the following expression for the pion-generating function Gt(^; z): 

G T (9 k ; z) = G (9 k ; z)(l - e^)2 2X ^R k \l + y k + R k f- 2 \ (40) 

where 



Rk = y/l-2x k y k + yl (41) 

z + (l-z) 2 sh 2 (9 k )ch 2 (9 k ) 
Xk z-{l-zfsh\9 k )ch\9 k y 

Vk = • 



. 0Ub z-(i-z)sh 2 {e k ) 



i + (i-z)*sh?(e k y 

and Go(9 k ; z) denotes the pion-generating function at the temperature T = 0: 

Go(9 k ;z) = [l + (l-z)sh 2 (9 k )]- 2X . (42) 

We observe that G is exactly the generating function of the NB distribution with a 
constant shape parameter 2A, and the average number of k-mode pions is equal to 

n(9 k ) = 2\sh 2 (9 k ). (43) 
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The vacuum value of the /c-mode thermal-density operator pr(0k) is used to obtain 
the /c-mode thermal eikonal function Qj<{9k) 

(0 k | Pr(0k) I fc ) = e- n *^ (44) 
= (l-e-^)G {9 k ;e-^). 

The total eikonal function is Qt(Y, B) = J2k ^r(9 k ). 

For the k-mode pion field in B space , we find the following average number 
and the dispersion: 

nr{9k) = n(Q k ) + n T k + jn(9 k )n T k, (45) 

2^ 3 4 

4(4) = d 2 Tk + d 2 {O k )[l + ^—n Tk + -n 2 Tk l 

where 

r2 -2 i - 

«T/e — n Tfc + n Tfc, 

d 2 (9 k ) = ^n 2 (9 k )+n(6 k ). (46) 

Two limiting cases are of interest, namely T — > and T — > oo. 
For T -> , we find 



^ 1 1 



(47) 



n 2 (^) 2A n(9 k y 
as it is to be expected from the NB distribution . However, the interpretation of 
this result is quite different. In our case, the parameter A is connected with the 
vacuum expectation value of the effective Hamiltonian H(Y, B). It has nothing to 
do with either the number of pion sources or the number of clans. Since SU(1, 1) is a 
dynamical symmetry group of our effective Hamiltonian, the parameter A also labels 
the positive discrete class of its unitary irreducible representations. It is important 

to observe that pions in the k mode are distributed according to the NB distribution 
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with a constant shape parameter 2A. The Wroblewski relation 

d(6 k ) = An(9 k ) + B (48) 

obtained from (47) gives the energy-independent coefficients A = (2A) -1 / 2 and B = 
(A/2) 1 / 2 . If A > 1/2, we have A < 1. 

The contribution from all the k modes in B space gives 

(49) 

— * 

where p(9 k ) = n(9 k )/n(Y : B). In this case, the coefficient A in the Wroblewski 
relation becomes energy and B dependent and is of the form 

k 

Since J2kP{@k) = 1 an d all p(9 k ) are positive functions of 9 k , the sum J2kP 2 (@k) is 

— * 

always smaller than unity. Therefore, A(Y,B) < 1 if A > 1/2. 

Finally, the summation over all impact parameters gives 

(^) 2 = /^Q(F,5)P 2 (F,5) + l)(^^) 2 -l] + ^. (51) 

This expression, when combined with our preceding analysis, suggests that the co- 
efficient A in the Wroblewski relation should be energy dependent and smaller than 
unity. 

For the temperature T going to infinity we obtain 

dU9 k ) ^ = 2 _ (1 + n^ y2 (52) 



n 2 Tk (9 k ) 

= l + th 2 (29 k ). 

This result shows that at very high temperature of the pion source the distribution 
of pions will tend to become chaotic if 9 k is very small. This will happen when 
the kinetic energies of the emitted pions are much larger than the corresponding 
coupling to the nucleon field, uj k g k {Y, B). 
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4 Conclusions 



In this paper we have proposed an intensity-dependent pion-nucleon coupling Hamil- 
tonian with 577(1, 1) dynamical symmetry. This Hamiltonian, within a multiparticle- 
production model of the AABS type in which the k - mode pion field is represented 
by the thermal-density operator, naturally explains the appearance of the NB mul- 
tiplicity distribution for pions in impact-parameter space. The shape parameter of 
the NB distribution is related to the vacuum expectation value of the Hamiltonian. 

The Wroblewski - type relation (1) has been obtained with the coefficient A 
that is energy dependent and smaller than unity if the vacuum expectation value of 
the Hamiltonian is larger than "zero-point energy" corresponding to A = 1/2. 

For T ^ 0, we have found a pion-generating function that may be used for 
obtaining all higher-order moments of the pion field. 

In our model, the k modes of the pion field are statistically independent and 
are therefore described by the factorized thermal-density operator. Correlations 
between different k modes are absent and, at this stage, our model cannot describe 
the emission of resonances. However, this can be remedied by adding a mode- 
mode interacting part to the Hamiltonian H(Y, B) [18]. We hope to treat this case 
elsewhere. 
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